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ABSTRACT 

Advances In flow measurement techniques In turbomachinery continue to be paced by the need to obtain 
detailed data for use In validating numerical predictions of the flowfleld and for use In the development 
of empirical models for those flow features which cannot be readily modelled numerically. The use of 
laser' anemometry In turbomachinery research has grown over the last 14 yr In response to these needs. 

Based on past applications and current developments, this paper reviews the key Issues which are Involved 
when considering the application of laser anemometry to the measurement of turbomachinery flowflelds. 
Aspects of laser fringe anemometer optical design which are applicable to turbomachinery research are 
briefly reviewed. Application problems which are common to both laser fringe anemometry (LFA) and laser 
transit anemometry (LTA) such as seed particle Injection, optical access to the flowfleld, and measurement 
of rotor rotational position are covered. The efficiency of various data acquisition schemes Is analyzed 
and Issues related to data Integrity and error estimation are addressed. Real-time data analysis tech- 
niques aimed at capturing flow physics In real time are discussed. Finally, data reduction and analysis 
techniques are discussed and Illustrated using examples taken from several LFA turbomachinery 
applications . 


INTRODUCTION 

The first application of laser anemometry to the measurement of turbomachinery flow fields was 
reported by Wlsler and Mosey In 1972 (Ref. 1). In the following 14 yr the quality and quantity of data 
generated by laser anemometer applications In turbomachinery has continued to Increase due to advances In 
optics, electronics, and computer hardware. This data has been used to Improve our understanding of the 
flow physical phenomena In turbomachinery and to validate numerical flow analysis schemes. 

Until recent years laser anemometer Investigations In turbomachinery have dealt with Isolated rotors 
since flow analysis techniques have until recently been confined to the steady, axlsymmetrlc flow regime 
which exists In the rotor relative reference frame for Isolated rotor configurations. However, several 
recent experiments have Involved the use of laser anemometry In studies of the periodically unsteady flow 
within the blade rows In single stage machines and between the blade rows In multistage machines. 

The application of laser fringe anemometry to measurements In turbomachinery environments Is reviewed 
In this paper. Example results are limited to axial-flow type compressor blading since most laser fringe 
anemometer experiments have Involved flow surveys In fans and compressors rather than turbines. Although 
radial-type turbomachinery has been surveyed using laser anemometry, the majority of these applications 
have Involved laser transit anemometry due to the generally superior ability of the LTA to make measure- 
ments In the narrow exit channels of centrifugal Impellers. A recent review of published results obtained 
In centrifugal compressors has been given by Kraln (Ref. 2). Additional LTA applications will be covered 

In detail In this symposium by Schodl and Elder (Refs. 3 and 4). 

FUNDAMENTALS OF LASER FRINGE ANEMOMETRY 

Basic Operating Principles 

The operating principles of laser fringe anemometers will be briefly reviewed below. In an LFA the 
laser output beam Is divided Into two equal power beams which are focussed to a common point In space 
which Is referred to as the measurement or probe volume. The crossing of the beams In the probe volume 
results In constructive and destructive Interference between the train of laser light waves contained In 
each beam. This Interference creates bright planes of light created by constructive Interference sepa- 
rated by dark planes caused by destructive Interference as shown In Fig. 1. The fringe planes are per- 
pendicular to the plane which contains the laser beams and are parallel to the beam bisector. As shown 
In Fig. 1 the spacing between bright fringe planes Is s = \/(2 sin K). A particle which crosses the 
probe volume scatters light at the fringe crossing frequency, f c . Note that the fringe crossing fre- 
quency Is determined solely by the component of particle velocity U x , which Is perpendicular to the 

fringe planes. The Uy and U z components carry the particle parallel to the fringe planes and do 
not therefore contribute to the fringe crossing frequency. Also note that while a rotation of the plane 
containing the laser beams about the beam bisector can be used to measure velocity components in direc- 
tions between the x- and y-dlrectlon, the line of sight velocity component U z cannot be measured. 

Figure 1 Illustrates three features which make laser anemometry (or LA) attractive compared to other 
velocity measurement techniques: 

(1) The system output, t.e., the doppler or fringe-crossing frequency, Is linearly related to the 

velocity 

(2) The technique Is not subject to drift since the proportionality between f c and U x is 

given by the laser light wavelength, x, and beam crossing angle, K, both of which are constant In 

time 


(3) The technique Is only sensitive to one component of velocity 
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Since each of the Incident laser beams Is circular In .cross-section, the actual probe volume shape 
Is an ellipsoid as shown In Fig. 2. The light Intensity distribution across each beam Is Gaussian In 
shape. The beam diameter, dg 2 , defined as the diameter at which the Intensity Is 1/e 2 of the peak 
Intensity at the center of the beam, Is typically used as a measure of the beam diameter. Using this 
definition of beam diameter, the probe volume diameter, d m , and length l m defined by the 1/e 2 Inten- 
sity level can be calculated as shown In Fig. 2. Note that the ratio of probe volume length-to-dlameter 
Is given by L/d = l/tan(K). Since the beam crossing angle K Is usually between 2 and 10°, the value 
of L/d Is on the order of 10 to 20. 

Another feature of Gaussian laser beams which Impacts laser anemometer optical performance Is that 
of beam divergence which Is Illustrated In Fig. 3. The two Important properties of the beam geometry 
shown In this figure are: 

(1) the laser light wave front radius R(z) Is Infinite only at the location of the beam waist, z = 0 

(2) the output beam waist diameter dg 2 . Is related to the Input beam waist diameter Dg 2 

through the formula dg 2 = 4\f/irD§ 2 . 

The first beam divergence property Implies that the probe volume fringes will not be parallel to one- 
another unless the Incident laser beams cross at the beam waist. If the fringe planes are not parallel, 
false levels of flow fluctuation will be Indicated since particles of equal velocity which pass through 
different parts of the probe volume will generate different fringe crossing frequencies due to the varia- 
tion In fringe spacing across the probe volume. A laser light collimator or mode matching lenses (see 
Ref. 5) can be used to Insure that the beam waist will be located at the probe volume location. 

The second beam divergence property shown In Fig. 3 Implies that the probe volume size can be con- 
trolled through control of the Incident laser beam diameter Dg 2 . Beam expanders can be used to 
Increase the beam diameter D§ 2 , which In turn reduces the beam diameter at the probe volume, d§ 2 . 

From the relations shown In Figs. 2 and 3 we see that expanding the beam diameter D§ 2 while holding 
the beam separation d constant results In a reduction In probe volume size and in the number of fringes 
In the probe volume. This results In an Increase In the power density In each fringe. This In turn will 
result In more photons being scattered from a given particle size or In the ability to generate the same 
number of scattered photons from smaller size particles. 

Optical Design Considerations 

The successful design of an LA optical system represents a trade-off between several conflicting 
factors. This Issue can best be addressed by considering an Illustrative example. 

The minimum variance which one may expect In a single velocity measurement has been shown In Ref. 6 

to be 
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c velocity of light 

h Planck's constant 

N number of fringes 

P laser power 

V velocity * 

W 0 probe volume radius 

n PMT quantum efficiency 

r 0 background light flux 

a s particle light-scattering cross-section 

Q solid angle of collection optics 

This relation Indicates that the measurement accuracy Is enhanced by: 

(1) reducing flare light and probe volume size 

(2) Increasing the collection optics solid angle, the laser power, and the number of fringes In the 

probe volume. 

Flare light can be reduced by using anti ref lectlve coatings on window surfaces and by optically mask- 
ing the collection optics. The solid angle of the collection optics can be Increased by using large aper- 
ture (l.e., small f-number) lenses. However, the lens diameter and price Increase for a given focal 

length as the f-number decreases and lens Imperfections become harder to correct. 

The probe volume size can be reduced by Increasing Dg 2 while the number of fringes In the 

probe volume can be Increased by Increasing the beam crossing angle K. However, one may reach a point 

at which there are not enough fringes In the probe volume or the fringe spacing Is too small. Referring 
to Fig. 2, consider the following example: d = 22 mm, Dg 2 = 1 mm, f = 200 mm. These parameters 

yield a probe volume diameter of d m = 131 ym containing 23 fringes with a fringe spacing of S = 4.68 ym. 

Commercially available counter-type signal processors can accurately measure fringe crossing frequencies 
up to 100 Mhz, but start to yield less accurate results for higher frequencies. The velocity which corre- 
sponds to a fringe crossing frequency of f c = 100 Mhz Is therefore a relevant parameter to consider and 

Is equal to V = f c .S = 468 M/sec for the present case. Increasing the Input beam diameter Dg 2 by a 

factor of 2 would reduce the probe volume diameter to 65 ym but would also reduce the number of fringes In 
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the probe volume to 14. This Is close to the eight fringes required by most counters and can lead to 
angle biassing errors If Bragg shifting Is not used. The number of fringes In the probe volume could be 
Increased back to the original value of 28 by doubling the Input beam spacing, d. However, this would 
change the fringe spacing to S = 2.35 ym and would decrease the velocity which corresponds, to a fringe 
crossing frequency of 100 Mhz to V = 235 M/sec. 

Measurement of Multiple Velocity Components 

There are many optical configurations which can be constructed for measuring one, two, and three 
velocity components simultaneously. Reference 7 contains a rather complete description of several one-, 
two-, and three-component configurations as well as a comparison of the relative merits of the various 
configurations. This Issue will be further addressed In this symposium by Boutler (Ref. 8) and will 
therefore be only briefly addressed below. 

In turbomachinery flowflelds the streamwlse and circumferential components of velocity V z and 
Vq respectively, are usually much larger In magnitude than the velocity component In the hub-to-shroud 
direction, V r . In addition, the radial velocity component generally lies along the optical axis of the LA 
system and therefore cannot be measured directly. Most LA applications to date have therefore been aimed 
at measuring the V z and V e velocity components. Increasing emphasis on secondary flow studies Is 
currently generating applications which require measurement of the V r velocity component as well. 

The simplest method of measuring velocity magnitude and flow angle Is to acquire measurements with a 
one-component LA system at two different fringe orientations. Use of a one-component LA system results 
In two limitations. First, the statistical error In velocity, V, and flow angle, a, measured with a 
single-component system Is greater than that resulting from a two-component ( system as will be described 
In the next section. Second, the magnitude of the turbulence components V z and V e can only be 
determined If measurements are taken In the z- and e-coordlnate directions or If measurements are made at 
three different fringe orientations. If measurements are acquired at two arbitrary fringe orientations, 
one can only determine upper bounds on V z and V 0 . Despite these limitations, single channel LA sys- 
tems have been used extensively In turbomachinery applications. The advantages of a single channel system 
are Its simplicity and the fact that the available laser power Is concentrated Into a single fringe 
system. 

Two-component LA systems can be Implemented by using two colors to create two separate measurement 
channels. Measurement channels can also be generated by using polarization separation or by using two 
Bragg cells to achieve frequency separation. If signal processor logic Is used to accept only those 
events for which a velocity measurement Is simultaneously made on both channels, the statistical error In 
calculating V and a is reduced relative to the error obtained when using a single-component LA system 
and the magnitude of the turbulence components V z , Vq, and V Z V 0 will be directly measured. 

Three-component LA systems can be Implemented by using three colors or two colors plus frequency sep- 
aration to create three measurement channels. These systems are quite complex. The accuracy with which 
such systems measure the radial velocity component Is directly related to the off-radial beam separation 
angle. These systems have transmitted beams which occupy a large solid angle and which therefore require 
relatively large windows for optical access to the flowfleld. In addition, the complex blade geometry 
found In most turbomachines may prevent such systems from having optical access to large areas of the 
flowfleld. 

It should be noted that one does not have to resort to a three-component LA system In order to obtain 
measurements of the radial velocity component. A technique In which the beams from a single channel LA 
system are deflected from the radial direction Is described In Ref. 9. This technique uses measurements 
obtained from two off-radial beam orientations to calculate the radial velocity component. As will be 
shown In the following section, the statistical error in the radial velocity component, V r , measured by 
this method will be greater than If V r were measured directly since this technique utilizes measure- 
ments made at two different times, 

A second approach to the measurement of radial velocities Is reported In Ref. 10. This approach 
utilizes the window configuration shown In Fig. 4 to enable direct direct measurement of the radial veloc- 
ity component when using the window labelled P'. While this approach does not require off-radial deflec- 
tion of the transmitted beams, It can only be used In regions which are upstream and downstream of the 
blade row. 

A third approach, which can be used to obtain radial velocity component measurements within a blade 
row, Is Illustrated In the right half of Fig. 5. The transmitted beams enter the flowpath at the blade 
stagger angle through an optical access hole located upstream of the blade row. Scattered light Is col- 
lected In an off-axis direction through a window located over the blade row. This optical arrangement 
provides measurements of V r and a combination of V z and V e . The separate V z and V 0 compo- 
nents are determined using the conventional approach shown on the left In Fig. 5. This approach Is appli- 
cable only to a single stage machine with no Inlet guide vane. Application of this approach to multistage 
machines would require the use of fiber optics to Introduce the transmitted laser beams Into the 
flowfleld. 

Statistical Measurements Errors 

In order to determine accurate estimates of the mean velocity and the turbulence properties of a 
flow, many Individual LA measurements must be acquired and averaged together. The mean and standard devi- 
ation estimated from the data are subject to both systematic and statlstlcal'uncertalnty. An analysis of 
these uncertainties Is given In Refs. 11 and 12. As shown In Ref. 11, the number of measurements required 
to establish a given level of confidence In the mean velocity Is given by 
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where V Is the true mean velocity, V m Is the measured mean velocity, (V'/V) Is the true turbulence 
Intensity, Z Is the confidence level, and N Is the number of measurements used to calculate V m . A 
value of Z = 1.97 corresponds to a 95 percent probability that the measured mean velocity lies within the 
range V(1 - C v ) < V m < V(1 + C v ). The expression for the confidence level for the rms velocity estimate 
Is 


V' 

m 


v 


V 


= c 


v' 



where V' Is the true rms velocity and V^ Is the measured value. Note that the confidence Intervals 
are Inversely proportional to V*- One therefore has to Increase the number of measurements by a factor of 
four In order to halve the confidence Interval. Also, note that C v Is proportional to the turbulence 
Intensity, V'/V, while Cy Is Independent of V‘/V. Since VW Is typically on the order of 0.1 or 
less, a given number of measurements always yields a much better estimate of the mean than of the rms 
velocity. These facts are Illustrated In Table I which shows the number of measurements required to 
establish a confidence level of 95 percent for various values of V'/V. 


An Important point addressed In Refs. 11 and 12 Is that the statistical error In calculating velocity 
components Increases when using uncorrelated data acquired at different fringe orientations. Data Is 
uncorrelated In a multichannel LA system If signal processor logic Is not used to enforce simultaneity of 
measurements from each channel. Uncorrelated data always occurs when using a single channel LA system to 
acquire data at multiple fringe orientations. The magnitude of the Increased statistical error can be 
Illustrated by considering how a single channel LA system Is used to determine orthogonal velocity 
components . 

Figure 6 Illustrates the calculation procedure by which velocity and flow angle are obtained with a 
single channel system. Measurement directions 1 and 2 shown In the figure are the directions normal to 
the probe volume fringes. Measured velocities V-; and V 2 , acquired at the known angles 6 ] and 63 , 
are used to solve for the unknown velocity magnitude, V, and flow angle a. Orthogonal velocity compo- 
nents V u and V v can then be determined from V and a. The analysis of Ref. 11 Indicates that the 
confidence levels In V u and V v are given by the expressions 



M u and M v are multipliers which Indicate the Increase In statistical error over the case In which V-| 
and V 2 are acquired simultaneously, l.e., M-| = M 2 = 1 when V-| and V 2 are acquired simultaneously. 

If we assume Isotropic turbulence, the expressions for M u and M v In terms of the geometry of Fig. 5 are 

2 2 i 1/2 

cos e 1 + cos ©2 

sin 2 (e 2 - e 1 ) 

Figure 7 Is a plot of these multipliers as a function of the mean fringe orientation angle and the magni- 
tude of the fringe angle difference, 02 - . These results Indicate, for example, that for the case of 

e-| = 10°, 62 * 70°, the statistical errors In determining V u and V v are Increased by factors of 1.217 
and 1.449, respectively when using data which Is not simultaneously acquired. The Impact of this 
Increased statistical error must be weighed against the Increased system complexity when deciding whether 
to use a single or multichannel LA system In a given application. 

Seed Particle Considerations 


M * 
u 


2 2 
sin 0 -| + sin 


sin ( 0 2 - 0 ^ 


l / c 


The tradeoffs between seeding the entire flowfleld or using a point source of seed are briefly dis- 
cussed below. In addition, the properties of a unique fluorescent seed material are described. General 
methods used to generate and size particles and considerations related to the particle size required to 
accurately follow the fluctuations In a given flowfleld will be addressed In detail In this symposium by 
Mel 1 1 ng (Ref. 13). 

When performing laser fringe anemometer measurements In a gas It Is usually necessary to Introduce 
seed particles Into the flowfleld for two reasons. First, the number of particles which are naturally 
present In atmospheric air In the 0.5 to 1.0 ym range Is generally not sufficient to yield adequate data 
rates. Second, by generating and Injecting seed particles Into the flow one can control the size distri- 
bution of the scattering particles and can therefore control to some degree the accuracy with which the 
seed particles follow the flow. 

When designing an LA experiment one must decide whether to seed the entire flowfleld (full coverage) 
or to seed only the stream tube which passes through the measurement volume (point Injection). Both 
methods are used In practice and the method chosen Is dependent on the particle generation rate of the 
seeder and on the flowfleld characteristics. As an example analysis of particle generation needs, let us 
consider an LA application which Involves flowfleld measurements In the turbomachinery environment sum- 
marized below: 


turbomachine diameter = 1 m 
hub/tip radius ratio = 0.7 
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through flow velocity = 200 m/sec . e ■' . 

volume airflow rate = 80 M 3 /sec 

The rate at which particles will cross the LA probe volume Is given by 

R = d ■ L* V*c 

where L and d are the probe volume length and diameter, c Is the particle concentration In the flow, 
and V Is the flow, velocity . 

Commercially available seeders are capable of maximum particle generation rates of 10 11 particles/ 
min. Coupling this rate with the volume flow rate of air through the machine yields a particle density 
of 2x1 0 7 particles/M 3 If full coverage seeding is employed. If typical probe volume dimensions of 
d = 100 ym, L = 1 mm are used, we arrive at a particle rate through the probe volume of 

R = d L V c * (10-* m)(10~ 3 m)(200 m/sec)(2xl0 7 /m 3 ) =• 400 particles/sec. 

The actual LA data rate will be less than this number since particles which cross the probe volume outer 
edges usually do not scatter enough light or cross enough fringes to yield a valid velocity measurement. 
The conclusion to be drawn from this particular example Is that point Injection of seed should probably 
be used. 

Point Injection of seed Is most easily accomplished through a tube placed upstream of the measurement 
volume as shown In Fig. 8. This tube must be placed far enough upstream to enable decay of the wake shed 
from the tube before the measurement point Is reached. In addition, the seed Injection tube Is usually 
carried In a radial and/or circumferential actuator to enable seed Injection of the stream tube which 
passes through the probe volume for arbitrary positions of the probe volume. 

A unique seed material which has specific advantages for turbomachinery LA applications Is the liquid 
fluorescing seed described In Ref. 14. This seed material Is formed by dissolving an organic dye, rhoda- 
mlne 6G, In a mixture of benzyl alcohol and ethylene glycol. When Irradiated with light from the blue or 
green lines of an argon-ion laser, seed particles composed of this dye solution fluoresce In the orange 
band of the spectrum, as shown by the absorption and emission spectra shown In Fig. 9(a). By placing a 

narrow band orange-pass filter In the LA collection optics In front of the PMT, one can optically stop 

light reflected from solid surfaces at the Incident light wavelength from entering the PMT. This enables 
one to make measurements near hub and endwall surfaces and within rotating blade rows. When not using 

fluorescent seed, the light reflected from solid surfaces, which Is orders of magnitude higher In Inten- 

sity than the light scattered by the seed particles, dominates the PMT signal as the surface Is 
approached, l.e., the signal to noise ratio of the PMT signal drops to zero. When fluorescent seed Is 
used, the signal to noise ratio remains constant as a solid surface Is approached. In addition, since the 
light from the seed particles Is emitted due to fluorescence rather than scattered, the Mle scattering 
phenomena does not apply and the intensity of light emitted by the particles Is Independent of direction. 
The penalty associated with the use of this technique Is that the light-emitting efficiency of the rhoda- 
mlne dye for a given size particle Is generally an order of magnitude lower than the backscatterlng effi- 
ciency of conventional (l.e., nonfluorescing) seed particles of the same size, as shown In Fig. 9(b). In 
addition, both blue and green wavelengths lie In the absorption band of rhodamlne as shown In Fig. 9(a). 
This technique therefore cannot be used with two-color LA systems. 

Alternate Laser Anemometer Methods 

In addition to laser fringe anemometry there are several alternate laser anemometry methods which can 
be used to obtain quantitative flowfleld measurements In turbomachinery. These methods are Fabry-Perot 
Interferometry; laser transit anemometry ( LTA) , (also referred to as the laser two focus (L2F) technique, 
the time-of-f light technique, and as the two-spot technique); and a hybrid method which combines features 
of both laser fringe and laser transit anemometers. 

The LFA technique measures the frequency difference between the light scattered from each Incident 
laser beam by a seed particle. In contrast, the Fabry-Perot technique directly measures the frequency of 
the scattered light by using a Fabry-Perot Interferometer as a scanning optical spectrum analyzer. Devel- 
opment of the Fabry-Perot technique for measurements of the llne^of-sight radial velocity components In a 
turbine cascade has been reported In Refs. 15 and 16 and will be covered In this symposium by Seasholtz 
(Ref. 17). 

The laser transit anemometer (LTA) will be covered In detail In this symposium by Schodl (Ref. 3). 

The principal advantages of the LTA technique over the LFA technique are: 

(1) the light Intensity In the probe volume Is higher than In a laser fringe anemometer which enables 
detection of smaller particles, 

(2) the sensitive length of the measurement volume In the llne-of-slght direction Is much shorter 
than In an LFA system which enables the acquisition of measurements near solid surfaces, 

(3) the LTA provides a high degree of sensitivity In flow angle measurements. 

The principle disadvantage of the LTA technique Is that It yields a low data rate due to the fact 
that the LTA measurement volume presents a much smaller "target" which a particle must hit In order to 
yield a velocity measurement (see Fig. 10). LTA data rates diminish rapidly as the turbulence level 
Increases In a flowfleld because the flow angle fluctuations generated by turbulence levels of 10 percent 
or more are large enough to prevent most particles from crossing both spots. 

A hybrid technique currently under development Is aimed at removing this limitation by expanding the 
circular LTA spots Into elliptic-shaped spots. The technique represents a trade-off between power density 
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In the probe volume, and probe volume target area. One approach to this technique Is described In Ref. 18. 
This technique will also be discussed In this symposium by Boutler (Ref. 8). 

APPLICATION OF LASER ANEMOMETRY TO MEASUREMENTS IN TURBOMACHINERY 

There are three major features of the turbomachinery environment which make the application of laser 
anemometry more difficult In turbomachinery experiments than In external aerodynamic experiments: 

(1) the close proximity of highly reflective surfaces (l.e., hub and shroud endwalls and blades), 

(2) flowpath geometry (curved endwalls and complex blade shapes), 

(3) the rotation of rotor blade rows. 

Each of these Issues will be addressed In this section. In addition, various data acquisition 
schemes will be compared and the advantages and limitations of each scheme will be discussed. Finally, 
the need to provide near real-time data reduction to support Intelligent execution of experiments will be 
discussed. 

Reflected Light Problems In Turbomachinery Applications 


Two separate problems related to reflected Incident light radiation arise In LFA turbomachinery 
applications. The first of these problems Is caused by reflections from window and hub surfaces when 
making measurements In the flowpath Inner and outer endwall regions. Reflections from the window surfaces 
can be reduced by using anti ref lection coatings on the window. Fluorescent seeding Is an effective means 
of eliminating the problem of poor slgnal-to-nolse ratio caused by reflection of the Incident laser beams 
from blade and endwall surfaces. 

While fluorescent seeding enables one to make measurements directly down to both the hub and shroud 
endwall surfaces, the proximity to a solid surface which one can achieve with nonf luorescent seeding 
before reflected light becomes a problem Is highly dependent on the type of surface In question and on the 
design of the optical system. It has been the author's experience that for optical systems with beam 
crossing angles on the order of 5° or less employing f/4 collection optics and no special stops In the 
collection optics, the closest approach to a hub endwall or a window In the shroud endwall Is on the order 

of 1 cm. For a machine with a diameter of 1 m and a hub/tip radius ratio of 0.8, a distance of 1 cm cor- 

responds to 10 percent of the blade span. 

* The second source of reflected light In turbomachinery applications is "blade flash" caused by pas- 
sage of the rotor blades through the Incident laser beams. The severity of the blade flash problem Is 
dependent mainly on the rotational speed of the blade row. The Incident laser beams "walk" up the surface 
of the blade as the blade rotates through them as shown In Fig. 11. As discussed In the next section, the 
blades appear to be "bent" In this figure due to spanwlse twist. The blade surface Is In the probe volume 
Itself for a very short time, near t * 0. However light can be reflected from the surface Into the col- 
lection optics at any point between time t = 0 when the probe volume first Intersects the blade surface 
and time t = when the beams are no longer blocked by the blade surface. Reflected light entering 
the collection optics during time 0 < t < tp can cause saturation of the photomultiplier tube (PMT). 

If this occurs the PMT requires a finite amount of time to recover and measurements which occur during the 
recovery time will be lost. For rotors operating at tip speeds on the order of 400 to 500 m/sec, blade 
passing frequencies. are on the order of 5 to 15 Khz, and time t^ Is on the order of 7 to 20 ysec. 

However for low speed machines this time can be at least an order of magnitude greater. It has been the 

author ' s ..experience that reflected light from blade surfaces does not cause PMT saturation In high speed 
machines but can cause problems In low speed applications. Fluorescent seeding therefore appears to be 
necessary In high speed applications only when measurements In the proximity of the hub or shroud endwalls 
are desired. 

Blade Geometry 

One feature of turbomachinery geometry which Impacts optical access to the flowfleld Is spanwlse 
twist of the blading as shown In Fig. 12. This twist leads to the "bent” appearance of hub-to-tlp blade 
sections such as AA and BB when viewed In a streamwlse direction. For blades which are radially stacked 
about the c.g. of each section this bending Is worst near the leading and trailing edges. If the Incident 
laser beams are constrained to enter the flowpath In the radial direction, shadowed regions occur as shown 
In which the tip blade section prevents optical access to the hub section. The extent of the shadowed 
region can be as great as 20 to 30 percent of the blade pitch at the hub. This shadowing can be elimi- 
nated If the LA optical system Is capable of orienting the Incident laser beams away from the radial 
direction as shown In Fig. 12. 

When acquiring measurements at design speed operating conditions, use of the blade design geometry 
Is generally an accurate method of calculating the Information required to relate the axial and circum- 
ferential measurement location to the blade surface location. However, at part speed conditions the blade 
geometry will vary from the design geometry due to deflections caused by reduced mechanical stresses. The 
principal deflection In axial-type blading Is an untwist or change In the blade stagger angle as shown 
schematically In Fig. 13(a). The amount of untwist varies from zero at the hub to a maximum at the tip. 
Axial blading also experiences some axial deflection due to thrust loading. In radial-type blading the 
primary deflection due to mechanical stress Is a "flowering" of the rotor disk, which causes a shift In 
the axial location near the rotor exit. These deflections can be accounted for by measuring the actual 
location of the rotor during data acquisition and correcting the design blade geometry as required. 

The blade geometry at the tip can be measured during data acquisition by positioning the probe volume 
at the tip radius. If seed Injection Is turned off, the fringes are oriented to measure the circumferen- 
tial velocity component, and the laser power Is reduced, then the LA signal processor will perceive the 
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blades passing through the probe volume as seed particles. Since velocity measurements will only occur 
over the blade tip, the blade circumferential location can be determined from the location at which mea- 
surements have occurred. The sensitivity of such a procedure to changes In the blade stagger angle can 
be controlled by the circumferential resolution of the LA measurements. Results obtained from an axial- 
flow fan using the above procedure are presented in Fig. 13(b). As described In Ref. 19 the measurements 
were acquired In a 51.3 cm diameter rotor using a resolution of 1000 measurement locations per blade 
pitch. The rotor has 22 blades which yields a circumferential resolution of 0.008 mm. The measured 
untwist was 1.96*0.04° compared to a finite element analysis prediction of 2.02°. These measurements also 
Indicated that the blade spacing around the wheel was uniform to within 2 percent. 

Window Design Considerations 

Optical access to wind tunnel models used In external aerodynamics experiments Is usually provided 
through flat windows mounted In the tunnel-walls. The use of flat windows In turbomachinery applications 
leads to distortions In the outer flowpath due to shroud curvature In both the circumferential and stream- 
wise directions, as shown In Fig. 14 for an axial-type blade. This flowpath distortion Is often greatest 
In the rotor tip region and can significantly alter the local tip clearance. This problem Is particularly 
acute In radial turbomachinery applications, where the shroud radius of curvature Is small In both the 
circumferential and the streamwlse direction. 

Flat windows are required for laser transit anemometer applications, since the highly focussed 
"spots" of a transit anemometer become highly defocussed If the window Is not flat. This defocussing 
reduces the power density In each spot to levels which are unuseable. Flat windows are not required for 
laser fringe anemometers. 

One solution to the need for flat windows In curved flowpaths Is the use of many small windows, as 
shown at stations 5 to 14 In Fig. 15. Such an approach Is the only one practical for radial-type turbo- 
machines If one desires to maintain proper tip clearances or to use laser transit anemometry. However, 
as seen In Fig. 15, the use of many small flat windows limits one's ability to acquire measurements at 
arbitrary streamwlse locations. 

For axial-flow configurations, In which streamwlse radii of curvatures are relatively large, curved 
windows can be easily fabricated by heating flat glass plates to the yield temperature of the glass and 
forming the glass over a mold. Such windows follow both the circumferential and streamwlse curvature of 
the flowpath and offer unlimited access to arbitrary axial and circumferential measurement points within 
the window boundary. 

Refraction at the air-glass Interface on window Inner and outer surfaces causes the following 
phenomena which Impact measurement accuracy: 

(1) The actual probe volume location Is translated relative to the location which would have occurred 
without refraction effects, 

(2) The Incident laser beams which generate the probe volume may uncross - 1.e. f the two beams will 
no longer lie In the same plane, 

(3) The beam crossing angle, which affects the calibration factor between fringe crossing frequency 
and seed particle velocity may change. 

Probe volume translation will always occur for both flat and curved windows, while changes In beam 
crossing angle only occur for curved windows. Beam uncrossing occurs when the solid angle between the 
Incident laser beam and the window surface Is not equal for each beam In a fringe anemometer. Beam 
uncrossing thus generally occurs for curved windows but can occur for flat windows when the beam bisector 
Is not perpendicular to the window surface. If the uncrossing distance at the probe volume location Is 
comparable to the beam diameter, the beams will no longer Intersect one another and the probe volume will 
be destroyed. A change In beam crossing angle can occur only across a curved window. The crossing angle 
changes when the change In surface normal orientations between the window outer and Inner surfaces are not 
the same for both laser beams. 

The magnitude of window refraction effects Is dependent on the beam crossing angle, the probe volume 
Immersion Into the flow, the window thickness and radius of curvature, and the orientation of the Incident 
laser beams relative to the window. These effects have been recently studied In detail for flat and 
cylindrical windows using an optical ray tracing technique which Is applicable to totally general window 
shapes (Ref. 20). The Impact of typical values of these parameters on the uncrossing problem Is shown In 
Table II. The conclusion which one draws from this example Is that relatively thin windows are required 
when curved windows are used. We generally use a window thickness of 3 mm in our research rigs, which 
have a diameter of 0.51 m. 

One very promising solution to the uncrossing problem Is presented In Ref. 21. The procedure 
Involves making an Interference hologram between a reference beam and a laser beam which passes through 
the window after passing through the desired measurement point as shown In Fig. 16. The hologram can then 
be used to Introduce an exact refraction-error correction Into an Incident laser beam which follows the 
path of the -reference beam used to create the hologram. The hologram deflects the beam onto the path 
required to pass through the window and the measurement point. The technique Is suitable for high laser 
power levels (on the order of 2 W) due to the high transmission efficiency of the hologram. The only 
disadvantage of this technique Is the need to make a hologram for each different location at which one 
wishes to pass through the window. 

Window Cleaning 

The use of seed Injection Into the flowfleld In order to enable LA measurements can result In the 
accumulation of seed particles on window Inner surfaces. This deposition Is due to large particles which 
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do not follow the flow and are therefore centrifuged out toward the shroud. In addition, turbomachinery 
component research rigs usually have lower quality oil seals than flight-certified engines, which can 
cause window contamination due to oil leakage even when seed material Is not used. An effective method 
of removing deposits from windows Is the Injection of a liquid solvent Into the flow through small holes 
located In the shroud upstream of the window. Window cleaning by this method Is especially effective In 
the region over the rotor tip, where the liquid Is scraped across the window by the rotating blades. The 

author has been able to wash windows In single-stage axial compressor applications for rotors with tip 

speeds on the order of 425 m/sec without reducing the rotor speed, which allows the operating point to 
remain unchanged during window washing. The use of liquids for on-line window washing may not be practi- 
cal for applications In higher pressure ratio centrifugal compressors and multistage axial compressors. 

In such machines the flowfleld temperature rise may be so high that window washing liquids would evaporate 
or cause high thermal stresses, when they coat a hot window. 

Measurement of Blade Rotational Position 

A laser anemometer measurement Is a random event triggered by the presence of a seed particle In the 
probe volume. When acquiring measurements In rotating blade rows one must therefore measure the blade row 
rotational position every time a velocity measurement occurs. There are four measurement methods which 
can be used to provide this rotational position: 

(1) Enable a measurement window once-per-rotor revolution during data acquisition. 

(2) Enable a measurement window once-per-blade passing during data acquisition. 

(3) Enable continuous measurement windows during data acquisition. 

(4) Record the time between measurements and assign measurement windows after data acquisition. 

Methods (1), (3), and (4) require a once-per-rev (OPR) timing signal from the rotor while method (2) 
requires a once-per-blade (OPB) timing signal. 

The generation of measurement windows for the first three methods Is shown schematically In Fig. 17. 
By definition a measurement window Is a time Interval used to classify velocity measurements according to 
the rotational position at which they occurred. All velocity measurements which occur within a given 
window are spatially averaged across the width of the window. It Is therefore advantageous to keep the 
measurement windows narrow. Since LA measurements are usually displayed graphically as discrete measure- 
ment points across the blade pitch. It Is Important to remember that each point represents a spatial 
average across the window width. 

Implementation of methods (2) and (3) for generating measurement windows requires the ability to sub- 
divide the time between once-per-blade (OPB) and once-per-rev (OPR) signals respectively. This require- 
ment Is met by generating pulse trains at frequencies which are Integral multiples of the OPB or OPR sig- 
nal frequencies. These pulses are accumulated In a counter which Is Initialized by each OPB or OPR 
signal. Whenever an LA velocity measurement occurs, the current counter contents are recorded with the LA 
measurement as a (velocity, clock count) data word pair. The clock count Is then used to determine the 
measurement window In which the velocity measurement occurred. 

Note that a variable delay after the once-per-rev signal Is needed In order to align the measurement 
windows with the blade passage. This effect Is Illustrated In Fig. 17 by delays Di and D 2 which 
are required for measurement axial locations ZMi and ZM 2 , respectively. Successful Implementation 
of methods (2) and (3) also requires the use of fairly high frequency pulse generators, as shown In the 
table below for two typical high speed compressor applications. 

FREQUENCY REQUIREMENTS FOR WINDOW-GENERATING 
PULSE TRAIN 

[NB - number of rotor blades; BPF - blade passing 
frequency; NWP - number of windows per blade 
passage; F - pulse train frequency; WW - window 
width In microseconds] 


Machine type 

rpm 

NB 

BPF, 

Khz 

NWP 

F, 

Khz 

WW 

Centrifugal 
Main blade 

36 000 

20 

12 

16 

192 

5.2 

Main and split- 
ter blade 


40 

24 

16 

384 

2.6 

Axial 

20 000 

40 

13.3 

50 

665 

1.5 


Use of constant frequency pulse generators will result In changes In the number of measurement win- 
dows between OPB or OPR signals due to rotor speed variations. One must therefore use a phase-locked loop 
or a tracking shaft angle encoder (Ref. 22) which Is capable of rapidly varying the pulse generation fre- 
quency In response to rotor speed variations. For high speed applications a digital shaft angle encoder 
Is preferred over a phase-locked loop, since the response lag in phase-locked loops Is generally too large 
to enable accurate tracking of rev-to-rev variations In rotor speed. 

The first method listed above Is the simplest, and Involves enabling a measurement window after a 
known time delay t has elapsed after occurrence of the OPR signal. A blade-to-blade velocity profile 
Can be achieved by collecting measurements for various delay times t. This method allows one to obtain 
measurements In different blade passages around the wheel. It Is the simplest of the four methods and was 
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therefore used In early applications of LA to turbomachinery (see Ref. 23 for example). While It Is sim- 
ple to Implement this method Is also Inefficient since the measurement window Is only opened once per 
rotor revolution. 


In method (2) a fixed number of measurement windows are generated between adjacent once-per-blade 
(OPB) signals. Since the OPB signals from each blade are all Identical, measurements which occur within 
a measurement window can come from any blade passage on the wheel. This method Is more efficient than 
method (1) since It opens a measurement window at all times. However, in order to measure the flowfleld 
In an Individual blade passage, the windows must be enabled for only the passage of Interest. In this 
case data can only be acquired during a time corresponding to (1/NB) of each revolution, where NB Is the 
number of blades. Reference 24 describes the use of this method In a laser transit anemometer system. 

Method (3) essentially opens measurement windows at all times and assigns a unique set of measurement 
windows to each blade passage. This technique therefore allows maximum data acquisition efficiency when 
one Is Interested In obtaining measurements In Individual blade passages. Use of this technique In a 
laser fringe anemometer system Is described In Ref. 22. 

The fourth data acquisition method Is similar to the method (3) In that It assigns a unique set of 
measurement windows to each blade passage. This method, which Is described In Ref. 25 for a two-bladed 
propeller, Is shown schematically In Fig. 18. During data acquisition, a trigger signal Is recorded once- 
per-rev along with the elapsed time between each LA measurement. The pltchwlse location of each measure- 
ment Is then determined after data acquisition Is complete by using the tlme-between-measurement data and 
the triggering pulse (or OPR) data. After the pltchwlse measurement location Is determined, the measure- 
ment is assigned to the appropriate window. Since measurement windows are defined by the user during 
post-run data processing, this method allows the user the freedom of changing the width and number of the 
measurement windows without rerunning the experiment. In contrast to methods (1) to (3), this method also 
enables the user to determine the circumferential location of each LA measurement within the measurement 
window. This Information can be used to assess the degree of spatial averaging which occurs when all 
measurements within a window are averaged together. 

The four approaches discussed above are by no means the only methods which can be used to determine 
the blade row rotational position for each LA measurement. An approach which Is a hybrid of the third and 
fourth approaches discussed above Is described In Ref. 26. This method Is similar to method (4) In that 
It determines the rotational position from the tlme-between-measurements and a triggering pulse. However, 
In this method the triggering pulses are generated by an encoder which generates up to 1000 pulses/rev 
(similar to method (3)), while In method (4) the triggering pulses are generated once-per-rev . 

Data Acquisition Methods 

The various methods of recording rotational position discussed above allow one to Implement several 
different data acquisition methods. The method which one chooses to use is dependent on the level of 
Information which one Is trying to extract from the measured flowfleld and on the level of flowfleld 
periodicity. 

Averaging Methods 

The simplest data acquisition method which enables one to determine the blade-to-blade distribution 
of velocity and flow angle within a rotating blade row can be referred to as a "blade average" mode of 
data acquisition. This data acquisition method can be Implemented using method (2) discussed above for 
recording the blade rotational position. When performing blade average data acquisition one generates 
measurement windows between successive once-per-blade signals. Measurements which are assigned to a given 
measurement window can originate from any blade passage around the rotor. When using this method one is 
therefore averaging together the flowfields which occur In each Individual passage, and all information 
concerning passage-to-pas sage flowfleld variations is lost. 


A second, more complex method of acquiring data can be referred to as a "window average" mode of data 
acquisition. This method can be Implemented using method (3) discussed above for recording the blade 
rotational position. All measurements which occur within a measurement window are known to originate from 
the particular blade passage In which the window lies. Window-average data acquisition therefore allows 
measurement of the blade-to-blade distribution of velocity and flow angle in Individual blade passages, 
which in turn allows one to assess the level of passage-to-passage flowfleld variations around the rotor. 
Although window average data acquisition provides more information than blade average methods, information 
concerning the distribution of the individual velocity measurements which occur within a measurement win- 
dow is still lost during data acquisition. Blade and window averaging methods average together all mea- 
surements which occur within a given measurement window. As shown In Ref. 9, If one records the number 
of measurements, N, which occur within each window, the sum of the measured velocities, EV^, 1 =1, N, 
and the sum of the squares of the measured velocities, I(V^) 2 , then one can calculate the mean and 
standard deviation of the measurements which occurred within the window using the following: 
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Note that only five words of information, N, IV, IV 2 , need to be recorded for each measurement window 
(one word for N which Is an Integer and two words each for IV and IV 2 which are real numbers). 
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In addition, the number of words stored Is Independent of the total number of measurements acquired 
during a run. 

As shown In Refs. 9 and 22, the window average mode of data acquisition can be an efficient method 
of obtaining good blade-to-blade spatial resolution while minimizing computer storage requirements. The 
LA system described In these references generates 50 measurement windows across 20 consecutive blade pas- 
sages. The system therefore records velocity measurements In 1000 measurement windows while using only 
5000 words of computer storage. As discussed In Ref. 22 and shown In Fig. 19, the window number (referred 
to as “shaft position'? In Ref. 22) Is used as an address In a stored data array. This scheme enables 
efficient real-time processing and graphical display of the measured data because the data Is assigned Its 
proper place In the data array as part of the measurement process. 

Data Capture Method 

The most complete method of data acquisition can be referred to as a "data capture" method. In the 
data capture method, each Individual velocity measurement which occurs within a window Is Individually 
recorded, and the velocities are not summed during data acquisition. This data acquisition mode therefore 
allows the user to examine the probability density distribution of the measurements which occurred within 
each window. 

There are several advantages and disadvantages to be considered when deciding which data acquisition 
mode to Implement. The data capture method lends Itself to direct memory access (DMA) transfer of data 
from LA counter processors to computer core memory since all measurements are recorded without any aver- 
aging or summing. However, this method can require large amounts of computer core memory and archival 
storage relative to the blade and window average methods since two words of Information (Integer repre- 
sentation of measured velocity, window number) must be recorded for each measurement. 

Consider for example a situation In which one wishes to record velocity measurements In 25 windows 
from blade-to-blade across 40 Individual blade passages. In addition, assume that one wishes to record 
100 measurements In each window and that the data Is uniformly distributed among the windows. Window 
average data acquisition In this case would require (25 windows/blade) by (40 blades) by (5 data words/ 
window), or 5000 words of storage. In contrast, data capture acquisition would require (25 windows) by 
(40 blades) by (100 meas . /window) by (2 words/meas . ) , or 200 000 words of storage. 

If DMA data transfer Is used, sorting of velocity measurements Into a data array which Is addressed 
according to window number cannot be done In real time, but must be performed either after data acquisi- 
tion Is complete or In parallel with data acquisition. Implementation of real-time data display when 
using data capture methods therefore requires use of a multitasking computer operating system or an Inter- 
ruption of data acquisition for calculation of real-time displays. These Issues are becoming less Impor- 
tant as the cost of computer hardware and software capability continues to fall. 

The ultimate acquisition approach would be to Implement data acquisition algorithms which can auto- 
matically screen the probability density distribution acquired In each window during capture mode acqui- 
sition. If the distribution Is well-behaved (l.e., nearly Gaussian), little Information Is lost. by simply 
recording the distribution mean and standard deviation. The output of such an Intelligent data capture 
method would therefore consist of probability density distributions for those measurement windows In which 
nonGausslan distributions were detected. 

Total Number of Measurements Required 

Two data acquisition Issues which remain to be addressed are: 

(1) how many measurement windows are required across a blade passage, 

(2) how many measurements are needed to establish statistical confidence In the data. 

The number of measurement windows required for accurate definition of the flowfleld Is dependent on 
the magnitude of the velocity gradients In the blade-to-blade direction. These gradients are averaged out 
across a measurement window during data acquisition since measurements which occur anywhere within a win- 
dow are effectively assigned to the center of the window. Figure 20 shows a blade-to-blade Mach number 
distribution measured near the tip of a transonic fan blade In which a two-shock passage shock system Is 
present. The measurements were acquired using 50 measurement windows across the passage. The Mach number 
distribution shown Is the average of the distributions which exist In each of the 17 blade passages In 
which measurements were acquired. In Fig. 20 the second shock Is clearly defined when N w ^ n( j 0w = 50. 

Also shown In Fig. 20 are the Mach number distributions which would result from using 25 and 12 windows 
across the blade pitch. Although the Impact of reducing N w -^ n( j ow to 25 Is minimal, a value of N w ^ n( j 0w 
of 12 leads to an almost total smearing of the second shock and does not show the Mach number Increase 
near the pressure surface. Note that method (4) discussed In the previous section concerning measurement 
of blade rotational position allows one to refine window resolution during post processing of the data. 

The other methods of measuring blade position require the user to fix the window resolution during the 
measurement process. 

As discussed In the section on statistical measurement errors above, the number of measurements 
required to establish a given level of statistical confidence In the mean velocity, V, Is related to the 
turbulence Intensity. The number of measurements required as a function of turbulence Intensity was sum- 
marized In Table I for various values of turbulence Intensity at the 95 percent confidence level. 

While a turbulence level of 5 percent Is not uncommon In the core flow regions of turbomachinery 
blade flowflelds, the turbulence Intensity level within wakes can rise as high as 20 percent. If one 
desires to establish a 1 percent confidence Interval In the mean velocity In a flow with a 5 percent tur- 
bulence level, then 100 measurements are required for each window. Acquisition of data across 1000 win- 
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dows would therefore require a^total number of 100 000 measurements to achieve the desired accuracy In 
each window If the number of measurements were uniformly distributed across all of the measurement win- 
dows. In reality the measurements are not uniformly distributed, as shown In Fig. 21 which shows the 
distribution of the number of measurements which occurred during the measurement of the Mach number dis- 
tribution previously shown In Fig. 20. One Is therefore left with one of two choices: 

(1) Collect data until the desired number of measurements are achieved In each window. 

(2) Collect the desired total number of measurements and accept an Increase In statistical error In 

those windows which have less than the desired number of measurements. 

For typical applications, In which the time for data acquisition Is of concern, the second choice Is 
more practical than the first. 

On-Line Data Monitoring 

The use of on-line data reduction and display can greatly Increase the effective use of an LA system 
in any application even If the data acquisition process Is not totally automated. On-line feedback from 
the measurement process can Indicate the presence of poor quality data and can also be used to ensure the 
capture of relevant flow features. 

Figure 22, adapted from Ref. 22, shows one possible real-time display. This display Is generated 
during a data acquisition process which uses 50 windows per blade passage across 20 blade passages. The 
display Is updated every 10 to 15 sec during a run and Is based on the data acquired up to that point In 
the run. The actual real-time display Is shown In the center of Fig. 22 and Is composed of the upper and 
lower plots shown In the figure. The upper plot Is a plot of the velocity profile averaged across a blade 
pitch, l.e., It Is calculated by averaging together measurements made In all of the Individual blade pas- 
sages. The blade suction surface Is on the left side of the plot and the blade Itself appears as the 
cross-hatched region on the right side of the plot. The lower plot In Fig. 22 Is a bar graph which shows 
the number of measurements acquired -In each of the 1000 measurement windows. The 20 gaps In the plot are 
caused by the 20 blades since no measurements are acquired as the blade passes through the Incident laser 
beams. The accumulation of dirt and fluorescent seed material on the surface of the ninth blade passage 
Is causing a large number of measurements due to light scattering from the blade surface. 

While a real-time data display can be used to assess the quality of the data at Individual survey 
locations, a rapid analysis of survey results Is needed to ensure the capture of relevant global flow 
features. Use of on-line data analysis can reveal Important data points which may have been missed during 
a survey. . These data points can then be acquired while the research rig Is still operating under the same 
conditions during which the survey was performed. An example of a case in which such a procedure would 
be useful Is shown In Fig. 23 for a test run In which the objective Is to map out the shock location In a 
transonic compressor rotor. The plot shown In the figure Is a streamwlse plot of relative Mach number 
created by plotting data acquired from several axial locations along line $-$ during a flowfleld survey. 
Results from the Initial survey, shown on the left, Indicate that the front passage shock lies somewhere 
between 10 and 20 percent chord and that the rear passage shock lies somewhere between 30 and 40 percent 
chord. If such a plot Is available during the research run, the LA operator can make additional surveys 
In these regions In order to more accurately determine the shock location. The results of such additional 
surveys, shown on the right In Fig. 23, Indicate that the actual front passage shock location Is at 
15 percent chord. 

Sophisticated data acquisition methods which feature automated data acquisition and on-line data 
analysis require relatively sophisticated minicomputers for successful Implementation. Use of 32-bit 
minicomputers with multitasking operating systems enables large data preparation, data acquisition, and 
data- reduction programs to be core-resident at the same time, enables sharing of data between such pro- 
grams, and enables concurrent execution. of such programs. Large minicomputers also provide the speed 
required to Implement on-line analysis of data acquired In a "data capture" mode. This analysis can pro- 
vide on-line Indication of flow features such as shock unsteadiness and wake vortex streets. Use of 
16-bit minicomputers requires that large data acquisition programs be subdivided and "overlayed" In memory 
due to core storage addressing limitations. The slower execution speed of these smaller computers also 
limits the researcher to the use of the more efficient averaging methods of data acquisition discussed 
above which are not capable of resolving certain flow features. Successful LA measurements can also be 
acquired In turbomachines using 8-blt microcomputers. The type of hardware and software required between 
these two extremes depends on the level of Information and the amount of Information which one wishes to 
extract from the LA measurements. 

ANALYSIS OF TURBOMACHINERY FL0WFIELDS USING LASER ANEMOMETER MEASUREMENTS 

A number of different data acquisition methods of varying levels of complexity were discussed above. 
It was shown that data acquisition schemes can be simplified by employing some level of data averaging 
during the measurement process. However, such averaging prevents the extraction of certain levels of 
Information from the flowfleld. In this section an analogous situation will be shown to exist for both 
data analysis and graphical data display methods. A variety of methods will be presented and Illustrated 
using examples from published results on laser anemometer applications to turbomachinery experiments. 

The Hierarchy of . Inf ormatlon Obtainable from LA Data 

■ The process of analyzing LA data can be schematically represented by the Informational pyramid shown 
In Fig. 24. The greatest amount of Information Is at the base of the pyramid. Successive application of 
averaging operators leads one to the peak of the pyramid. However, detailed Information on flow features 
Is lost as one moves to higher levels of the pyramid. 

When performing; an analysis of the variance In LA measurements acquired from flowflelds generated 
within rotating turbomachinery one must be mindful of the periodicity present In the flowfleld. In sta- 
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tlonary, ergodlc flows, one usually equates the level of measurement variance with the level of random 
flow fluctuations due to turbulence. As will be shown below, there Is at least one contribution to the 
measurement variance which Is not related to turbulence at each level of the pyramid. 

The principles discussed above will be Illustrated by analyzing an actual LA data set. This data set 
was acquired one-half chord downstream of an axial flow fan with a one-component LA system using 50 mea- 
surement windows between each blade. Measurements were recorded In 17 of the 22 fan blade passages. The 
total number of measurements recorded was 30 000, which yielded an average of about 40 measurements In 
each window. The velocity component measured Is that which Is In the streamwlse direction outside of the 
blade wakes. 

Analysis of Measurement Mean and Variance 

The base of the Informational pyramid contains the Individual velocity measurements acquired In a 
capture mode of data acquisition. This data can be visualized as a two dimensional array V(1,j). The 
measurement window number, 1, ranges between 1 and NW = NWP*NB, where NW Is the total number of mea- 
surement windows, NWP Is the number of windows per blade passage, and NB Is the number of blade pas- 

sages surveyed. In the example data set being used here, NW = 50, NB = 17 and NWP = 850. The Index J 
ranges between 1 and N(1), the total number of measurements In window number 1. As shown In Fig. 25, 
this data can be presented In the form of a velocity probability density distribution (p.d.d.) for each 
measurement window by classifying the N(1) velocity measurements Into velocity bins of fixed width. 
Although p.d.d. ' s can be generated for each measurement window, there are not enough measurements In an 
Individual window In this data set to given an accurate description of the blmodal character of the 
p.d.d. 's which lie In the blade wake (point A). The p.d.d. 's shown In Fig. 25 therefore contain measure- 
ments acquired within the measurement windows located at points A and B In all 17 measured blade passages, 

and are therefore "blade average" p.d.d.'s. As shown In Ref. 27, the blmodal character of the p.d.d. at 
point A Indicates that a vortex street Is present In the blade wake. 

Time averaging of the data acquired In each window leads to the second level of the Informational 
pyram1d,_wb1ch contains the time-average circumferential velocity distribution between Individual rotor 
blades, V(1). This velocity distribution can be directly acquired using a "window average" data acquisi- 
tion technique or can be calculated from capture mode data using the averaging operation 

V(1) = j^yy £ V(1,J) 1 = 1 NW (1) 

Note that V(1) Is the statistical mean of the probability density distribution of the measurements from 
window number 1 . 

The time-averaging operation also allows one to calculate the standard deviation, o r (1), of the 
probability density distribution In window number 1. This standard deviation Is the lowest level of rms 
velocity Information available (see Fig. 24) and can be calculated using the formula 



Once the data Is time-averaged, detailed Information concerning the p.d.d. of the Individual measure- 
ments which occurred In the measurement window Is lost. However, the general level of broadening In the 
p.d.d. Is Indicated by the magnitude of cr r ( 1 ) . In addition. Information on passage-to-passage varia- 
tion In the flowfleld Is still present In the data. Examples of the circumferential velocity distribution 
given by Equation (1) are shown In Fig. 26 across 2 of the 17 measured passages. The flags at selected 
points Indicate the magnitude of the standard deviation o r . Note that the standard deviation rises In 
the wake as one would expect due to viscous mixing. However, the standard deviation alone does not Indi- 
cate the vortical nature of the wake flow as Inferred by the p.d.d. shown In Fig. 25. This Information, 
which was contained In the p.d.d.'s has been lost due to averaging. 

Now let us consider the factors which contribute to the variance a r . In a steady, laminar flow, 
all measurements In a probability density distribution would fall into the same velocity bln. In reality, 
the p.d.d. Is broadened by the following mechanisms: 

(1) Random turbulent fluctuations In the flow, 

(2) Flow unsteadiness which occurs at frequencies which are not Integral multiples of the rotor 
rotational frequency (such as rotor shock oscillations and trailing edge vortex shedding), 

(3) Averaging of velocity gradients across the width of the measurement window, 

(4) Flow variations caused by rotor speed drift during data acquisition, 

(5) Individual LA measurement errors. 

The contributions of Items (2) to (5) must be subtracted from the standard deviation cr r (1) cal- 
culated from the p.d.d. In order to obtain the level of turbulent flow fluctuations. The LA user can 
control or can estimate the magnitude of some of the elements In Items (2) to (5) above. 

The contribution of Item (3) can be minimized by reducing the physical size of the measurement win- 
dow. The contribution of Item (5) can be controlled by selecting the fringe orientation angles and the 
number of measurements so as to minimize the statistical measurement error. Individual LA measurement 
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errors arising from particle tracking errors can be minimized by generating uniformly small seed parti- 
cles. As discussed in Ref. 9, the combined contributions of Items (3) and (4) are on the order of 1 to 2 
percent for a practical application. 

Application of a spatial ensemble averaging process to the circumferential velocity distribution 
V(1) leads to^the third level of the Informational pyramid, which contains the ensemble-averaged velocity 
distribution 7(k).- The ensemble averaging operation Is given by 

NB-1 

V(k) = jjjg X V(k +■ m.NWP) k=l , . . . , NWP (3) 


This ensemble averaged velocity distribution can be directly acquired using a "blade average" data 
acqul si tlon' .technique In which NWP windows are generated between successive blade passing signals. The 
V(k) velocity distribution contains Information on the circumferential variation In the flowfleld across 
an "average" blade passage. Detailed Information on the variation of the flowfleld between Individual 
blade passages Is lost. This Information Is lumped Into the standard deviation a g (k) at each point 
In the ensemble-averaged velocity distribution, given by 


cAk) = 


NB 


NB-1 

£ 

m=0 


[V(k) - V(k + m. NWP ) y k=l, 


, NWP 


(4) 


Fig. 27 Is an Illustration of the ensemble-averaged velocity distribution calculated from the example 
data set by averaging across the 17 measured blade passages. The flags shown on selected points Indicate 
the magnitude of the rms velocity o g . 

This rms velocity arises from passage-to-passage flow variations. This geometric contribution to rms 
velocity fluctuations can be seen In Fig. 28, where the Individual velocity distributions across the 17 
measured blade passages are plotted to the same scale. 

The fourth level of the Informational pyramid contains the pi tchwl se-averaged velocity . 7. This 
velocity can be obtained from the spatially ensemble-averaged velocity using the averaging operation 

, NWP 

5 * SSP £ v < k > < 5) 

The standard deviation of V Is calculated using 


°p = ~ ? < k >] < 6 > 


The values of 7 and o p calculated from the ensemble-averaged velocity distribution are Included 
In Fig. 27. This Information can be directly acquired by a simple LA system which contains no capability 
for measuring rotational position of moving blade rows. All details concerning the blade-to-blade flow- 
field are lost at this level of averaging. Data of this type Is of Interest when comparing LA measure- 
ments to those obtained with aerodynamic Instrumentation such as total and static pressure probes and 
thermocouples whose frequency response is well below the blade passing frequency. The time-average 
velocity at the measurement point Is also used when measuring the steady state flowfleld In stationary 
blade rows. 


At the top of the Informational pyramid, the variance In the pitchwl se-average velocity, ©t, con ~ 
tains contributions from all lower levels of the pyramid. Several authors (Refs. 28 to 30) have addressed 
this Issue In the past using a model similar to that shown In Fig. 29. These Investigators have shown 
that the total rms velocity at the measurement point, a is related to the random fluctuation level and 
the periodic fluctuation level by the relation 


2 

°t 


0 



0 


2 

P 


These Investigators have assumed that the flowfleld Is the same In each blade passage. If however, 
there are flowfleld variations from passage-to-passage, such that the velocity distribution 7(1) in each 
passage is not Identical to the ensemble-averaged distribution ?( k), an additional contribution oq 
must be Included: 

2 2 2 2 
a t = °r + °p * °g 

It should be noted that calculation of the ensemble-averaged profile, V^k), can be performed directly 
from the 1 ndl vldual _veloc1 ty measurement array V(1J). It should also be noted that the pltchwlse- 
averaged velocity, 7, Is Identical to the arithmetic average of all measurements In the V(1,j) array. 

The calculation of V(k) In terms of V(1) shown In equation (3) and the calculation of 7 In terms of 
\f( k) shown In equation (5) was done In order to more clearly show how Information Is lost during each 
averaging operation. The following formulas can be used for direct calculation of 7(k) and 7: 


f NB-1 N(k+m-NWP) 1 

£ 

T V(k + m • NWP, j ) 

Lm=0 

3=1 _ _ 1 


~NB-1 


£ N ( k + m.NWP) 


_m=0 


k=l NWP 
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Care must be taken when computing all of the averages discussed above, since the velocity Is a vector 
which has both magnitude and direction. This point can be seen more clearly by considering the averaging 
of the two vectors shown In the diagram below: 



The actual average velocity obtained by vector addition Is of course zero. However, If one averages 
the velocity magnitude and angle separately one gets 


<vv 


(Bj ♦ 0 2 ) 


= 100 

(45 - 135) 


= -45 


which Is incorrect. 

The proper method of calculating average values of any velocity triangle parameter Is to average the 
primary velocity components V x and V y first, and then calculate the desired quantity. In the 
above Illustration this would lead to 

(U V (v xl * V x2_l (71.7 - 71.7) 

1 Vav = 2 2 ~ u 


(V, 


< V v1 * V v2* _ (71.7 - 71.71 


y av 


= 0 


The proper method of calculating the ensemble-averaged flow angle, B(k), for example Is therefore: 


N B- 1 

V (k) * X V (k+m.NWP) 
1 m=0 1 


1 NB-1 

V 0 (k) • -g Z V e (k+m.NWP) k=l , . . . , NWP 


B(k) 


= ARCTAN 


Vk) 

v z (k) 


k=l NWP 


This same procedure applies to other velocity triangle parameters such as relative and absolute 
velocity magnitude and relative flow angle. 

The preceding discussion Indicates that the level of Information which one wishes to extract from the 
flowf jeld can. be used to determine both the type of data acquisition method used and the level of averag- 
ing employed during data post-processing. Although Information Is lost through averaging as one proceeds 
up the pyramid shown In Fig. 24, substantial savings In data acquisition time, data processing time, and 
data storage space can be realized at the upper pyramid levels. For example, let NR be the number of 
measurements required to describe the. velocity to a certain level of confidence at each point In the cir- 
cumferential velocity distribution V(1). In keeping with the example used above, let the total number 
of windows be NW=850, corresponding to 50 windows per blade across 17 blades. The total number of meas- 
urements required Is therefore NT=NR*NW = 85Q*NR. If, however, one Is only Interested In the ensemble 
averaged velocity \f or the time average velocity 9, then the same level of confidence requires only 
50*NR and NR measurements respectively. 


Continuing with the above example, one cannot hope to assimilate all of the Information contained In 
the p.d.d.'s_for each of the 850 measurement windows. Consideration of the circumferential velocity dis- 
tribution, 7(1) across 17 blade passages at each axial, radial survey location within the blade row can 
also be quite tedious. Therefore, In practice one usually resorts to detailed consideration of one of the 
following : 

(1) The distribution of 7(1) across a few Individual blade passages If passage-to-passage flow 
variations are large. 

(2) The distribution of \T(k) across the ensemble-averaged blade passage If passage-to-passage flow 
variations are small. 


Graphical Methods of Oata Presentation 

Several different methods of presenting the data accumulated along the circumferential measurement 
line swept through a rotating blade row at one axial, radial survey location have been presented In the 
preceding sections. Additional graphical methods are required when attempting to visualize global flow 
features . 
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Shock wave locations In transonic blade rows can be determined from a series of blade-to-blade plots 
of data acquired at different streamwlse locations. An Independent determination of shock location can 
also be achieved by streamwlse cross-plotting of data obtained at a constant pltchwlse distance. These 
combined procedures are Illustrated In Fig. 30. Although the seed particle velocity may lag the gas 
velocity for a short distance downstream of the shock, the point at which the measured velocity first 
begins to decrease provides a consistent measure of the shock front location. 

In contrast to the one-dimensional slices through the flowfleld offered by blade-to-blade and stream- 
wise plots, contour plots offer the user a two-dimensional view of the flowfleld. An example of a contour 
plot of the Mach number distribution In a transonic axial-flow fan Is shown In Fig. 31. Contour plots are 
often used when comparing LA measurements to results generated by numerical flow analysis codes. 

-Changes In the magnitude of a flow parameter are more easily discerned from one-dimensional plots 
such as those In Fig. 30 than from contour plots. This Is due to the fact that one must visually Inte- 
grate across contour lines to determine magnitude changes when viewing contour plots. As a result, 
contour plots tend to be more qualitative than quantitative. This situation can be Improved by three- 
dimensional plotting which adds parameter magnitude to contour plots as a third dimension. The use of 
varying viewing angles, perspective, and hidden line techniques can greatly enhance the Information 
obtainable from three-dimensional plots, as shown by the examples In Fig. 32. The valley-like features 
ahead of the blade row represent the Mach number variation caused by bow waves from adjacent blades. The 
passage shock location Is clearly shown In the lower part of Fig. 32. 

Example Results from LA Applications In Turbomachinery 

Published results from LA Investigations will be briefly reviewed below In order to demonstrate the 
capabilities of the LA technique. The examples will be limited to axial-type turbomachines. A more com- 
plete list of examples from axial-flow turbomachinery experiments Is given In Ref. 31. An excellent sum- 
mary of LA applications In radial-type machines has recently been given by Kraln (Ref. 2). In addition, 
results obtained In radial-flow turbomachinery will be discussed In this symposium by Schodl and Elder 
( Refs . 3 and 4) . 

Passage-To-Passage Flow Variations 

An Ideal turbomachine would produce an Identical flowfleld In each blade passage In the blade row. 
However, manufacturing tolerances limit the degree to which real machines approach this Ideal goal. 
Passage-to-passage flow variations measured In two transonic fan rotors are shown In Fig. 33 In the form 
of blade-to-blade distributions of relative Mach number In Individual blade passages. The blade suction 
and pressure surfaces are denoted by SS and PS, respectively. The results shown In Fig. 33(a) were 
obtained In a 550 m/sec tip speed fan with an aspect ratio of 2.87 which had a part-span damper. The 
results shown in Fig. 33(b) were obtained in an undampered 1.56 aspect ratio fan with a tip speed of 
429 m/sec. The Mach number distributions shown In Fig. 33(a) for blade passages 1 and 15 encompass the 
range of those measured In each of the 38 Individual blade passages. Note that the flowfleld In passage 
15 Is representative of an oblique shock while the flowfleld In passage 1 Indicates the presence of a 
normal shock. Analysis of data from additional fans has led to the conclusion that the large variations 
shown In the dampered fan flowflelds arise due to variations In blade geometry caused by manufacturing 
variations In the dampers. Another feature which appears in Fig. 33 Is the Improved flowfleld uniformity 
at the near stall operating condition shown In Fig. 33(b). This phenomena has also been observed In data 
from other LA experiments and appears to be a general behavior of transonic fans. Large passage-to- 
passage flowfleld variations have two Important Impacts: 

(1) They generate an Increased level of apparent turbulence to downstream blade rows as represented 
by the rms velocity o g discussed In the previous section. 

('2) They make numerical flowfleld analysis difficult because there Is no single blade passage In 
which the flow can be considered to be representative of the flow throughout the blade row. 

Unsteady Flow Features In Isolated Rotors 

The capture mode of data acquisition allows detection of flow fluctuations with frequencies which are 
not integral orders of the once-per-rev frequency. One example of this capability Is the detection of 
vortices In the blade wake as Inferred f-rom the double-peaked probability density distribution shown In 
Fig. 25. Another example of this capability is shown In Fig. 34. which Is taken from Ref. 19. The 
blmodal nature of the probability density functions shown In Fig. 34 Indicates an oscillation of the rotor 
passage shock about the mean location shown In the upper part of the figure. The blmodal p.d.d.'s 
obtained In windows C and 0 are due to the fact that these windows lie upstream of the shock part of the 
time and downstream of the shock part of the time. 

Blade Row Interactions 

The unsteady flow features just discussed occurred In Isolated rotors where the flowfleld should 
nominally be axtsymmetrlc and steady In the rotor relative frame of reference. When a rotor Is followed 
by a stator, the flowfleld Is unsteady and nonaxl symmetric In both the relative and absolute frames of 
reference. An Investigation of the periodically unsteady flow through a stationary blade row therefore 
requires a circumferential survey across one blade pitch In the stationary blade row. The array of survey 
locations used during the experiment reported In Ref. 32 to Investigate the flow through a stator located 
downstream of a compressor rotor Is shown In Fig. 35. Note that a total of 128 survey points were used 
In order to survey the flow at 15 axial locations, which represents an eightfold Increase In the number 
of survey locations which would have been required If the flow was axl symmetric . Data acquired at each 
pltchwlse location must still be tagged with the rotor rotational position. By plotting the velocities 
measured at each axial and circumferential location for a given rotor rotational position one obtains a 
snapshot of the periodically unsteady flow through the stator. Figure 36 Illustrates this procedure for 
three different pltchwlse locations of the rotor relative to the stator. Turbulence Intensity Is chosen 

/ 
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as the plotted parameter In this case In order to clearly define the rotor wake location. The rotor wake 
Is sheared as It convects through the stator passage due to the fact that the velocity Is higher on the 
suction side of the passage than on the pressure side of the passage. A similar analysis of turbine vane 
wakes passing through a turbine rotor has been reported In Refs. 30 and 33. 

The above Investigations were performed In single stage machines. Measurements acquired within a 
multistage axial flow compressor using the laser transit technique have been reported In Ref. 34. The 
application of the LFA technique to a high speed multistage axial compressor will be discussed In this 
Symposium by Williams (Ref. 35). 


CONCLUDING REMARKS 

During the last 14 years laser anemometry has been shown to be a viable technique for obtaining 
detailed Internal flow measurements from within the hostile turbomachinery environment. Advances In 
optics, electronics, and computer technology during this time have led to the availability of commercial, 
off-the-shelf laser fringe and laser-transit anemometer systems for use In turbomachinery Investigations. 
Future advancements In LA applications will therefore depend on Improvements In the following areas: 

(1) The efficient acquisition, analysis, and display of the large amount of data which can be pro- 
duced during LA applications, 

(2) Measurement of the radial velocity component which usually lies In the llne-of-slght direction 
of the LA system, 

(3) Acquisition and analysis of data from within multistage turbomachines, 

(4) Use of flow visualization techniques such as holography and trace gas Injection to scope global 
flow features In order to reduce the number of physical locations at which LA measurements must be 
acquired, 

(5) Accurate measurement of turbulent and periodically unsteady flow fluctuations In turbomachinery 
In order to accurately assess the relative Importance of random and coherent fluctuations, 

(6) development of optical techniques for measurement of thermodynamic properties which can be used 
In conjunction with the LA technique to yield simultaneous measurement of velocity and density, pressure, 
or temperature. 
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TABLE I. - NUMBER OF MEASURE- 
MENTS REQUIRED FOR 95 PERCENT 
CONFIDENCE LEVEL IN MEASURED 
MEAN VELOCITY AND TURBULENCE 
INTENSITY 

[N \j = number of measurements 
required for 95 percent con- 
fidence that the measured 
mean velocity is within Cy 
percent of the true mean 
velocity; N v , = number of 
of measurements required for 
95 percent confidence that 
the measured turbulence 
intensity is within C y < 
percent of the turbulence 
intensity; V'/V = turbu- 
lence intensity.] 


Mean Velocity Estimate 

Desired 

Turbulence level. 

accuracy. 

V'/V, 

Dercent 

c V ' » 



percent 

1 5 

10 20 

0.5 

16 400 

1600 6400 

1.0 

4 100 

400 1600 

2.0 

1 25 

100 400 

5.0 

4 

16 64 

10.0 

1 

4 16 

Turbulence Intensity 



Estimate 


Desired Number of 

accuracy, measure- 

C V 

, ments, j 

J percent | 1 

Ny » 

O.S 

80 

000 

1.0 20 

000 

2.0 5 

000 

5.0 

800 

10. C 


200 




TABLE II. - LASER FRINGE ANEMOMETER BEAM REFRACTION EFFECTS ACROSS 
A CYLINDRICAL WINDOW 


W i ndow 
thickness, 
in. 

Beam 

orientation 
angle, 6 
degrees 

Beam crossing 
angle at the 
probe volume, 
degrees 

Radial shift 
in p.v. loca- 
tion, aR 
in. 

Beam uncrossing 

distance 

Microns 


Percent 
of probe 
volume 
diameter 


probe volume radius = 

9.5 in. 



0.125 

0 


1 

pi 

m 

0 


45 


1 



2.4 

.250 

0 



o 


0 


45 

1.405 

.080 

6.63 


5.3 


probe volume radius = 

5.0 in. 




0.125 

0 

1.4000 

0.043 

0 

0 

0 


45 

1.4015 

.027 

19.7 

1.87 

15.8 

.250 

0 

1.4000 

.086 

0 

0 

0 


45 

1.4025 

.053 

39.6 

3.76 

31.7 
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Figure 1. - Fringe model of a laser anemometer. 
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Figure 2. - Geometry of a laser fringe anemometer. 
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Figure 3. - Properties of a Gaussian beam. 
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Figure 4. - Measurement of radial velocity components downstream of a 
compressor blade row. yi, V0 measured using window P. VR mea- 
sured using window p'. 







(a) Conventional system. 



ROTATION 


(b) Radial flow system. 

Figure 5. - Measurement of radial velocity components 
within a blade row using off-axis light collection. 
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Figure 7. - Increase in statistical error due to use of 
uncorrelated measurements in calculating velocity 
components. 



Figure 8. - Point injection of seed into streamtube 
passing through the measurement point. 
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(a) Emission and absorption spectra. 
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(b) Comparison between fluorescence emission and light 
scattering versus particle size. 

Figure 9. - Characteristics of rhodamine 6G fluorescing dye. 
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Figure 10. - Comparison of typical laser-fringe and laser-transit anemometer 
probe volumes. 




Figure 11. - Incident light reflection from blade surfaces. 



(a) Axial view. (b) View looking downstream. 

Figure 12. - Existence of shadowed regions caused by blade twist. 






Figure 13. - Comparison of measured and predicted 
untwist at the tip of an axial flow fan at a tip speed 
of 429 m/sec. 



(a) Tip clearance distortion due to casing 
curvature. 



(b) Tip clearance distortion due to flow- 
path convergence. 

Figure 14. - Flowpath distortions caused by flat 
windows. 



Figure 15. - Use of small flat windows in radial tur- 
bomachines to minimize flowpath distortions. 
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(b) Holographic refraction correction. Incident laser 
beams propagate through the hologram in the 
direction taken by the reference beam during holo- 
gram construction and cross at the desired probe 
volume location. 


Figure 16. - Holographic correction of window refraction effects. 













Z MEAS 


(Z M>! 


(ZM) Z 



OPR 


(c) Generation of continuous windows. 


Figure 17. - Three methods of generating measurement windows for LA applications in rotating blade rows. 
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(a) Exact circumferential location of each measurement calculated 
using elapsed-time-between-measurements data and trigger 
pulse data. 



(b) Measurements assigned to proper window during post-run 
processing. 

Figure 18. - Assignment of measurement windows during post-test 
data reduction. 




V(X) = ONE VELOCITY MEASUREMENT AT SHAFT POSITION X 
SP « SHAFT POSITION 


Figure 19. - Data storage scheme used during "window average" mode data acquisition. 
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Figure 20. - Effect of window width on resolution of 
blade-to-blade velocity gradients. 
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Figure 21. - Example of nonuniform distribution of the 
number of measurements which occur in each mea- 
surement window. 







Figure 22. - Real-time display used in the NASA-Lewis 
system during compressor research testing. 




(a) Streamwise distribution (view S-S) of relative 
Mach number. Initial survey - every 10 % 
chord over first 40% chord. 



(b) Streamwise distribution (view S-S) of relative 
Mach number. Refinement of shock location 
using additional survey locations. 

Figure 23. - Use of on-line data analysis to capture 
shock location in a transonic compressor. 
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Figure 24. - Schematic representation of information available from LA data. 
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(a) Velocity distribution measured across 6 of the 17 measured 



blade passages. Average velocity = 678. 2. 



(c) p.d.d. for data acquired at point B in all 17 measured 
blade passages. Average velocity = 773. 0. 

Figure 25. - Blade-averaged probability density distributions. 
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- Figure 28. - Contribution of passage-to-passage flow 
variations to the total rms velocity level. Passage- 
to-passage flow field variations generate rms 
velocity Og. 
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Figure 29. - Model used to illustrate various contributions to the 
total rms velocity level. 
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(a) Blade-to-blade distribution (view B-B) of relative 
Mach number. 



(b) Streamwise distribution (view S-S) of relative Mach 
number. 


Figure 30. - Use of pitchwise and streamwise Mach number 
distributions to determine shockwave locations. 
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Figure 31. - Contour plot of relative Mach number distribution 
in a transonic fan. 
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(a) Dampered fan, 550 m/sec. 
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(b) Undampered fan, 429 m/sec. 

Figure 33. - Passage-to-passage flow variations in two transonic 
axial flow fans. 
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(b) Velocity probability at points near the passage shock. 

Figure 34. - Analysis of rotor passage shock oscillations using velocity probability density 
distribution. 





\ 


ROTOR 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


o o o 
0^000 
o ooo o o o 
o ooo ooo 
o ooo ooo 
o ooo ooo 
6 ooo ooo 
o ooo 
o 
o 
o 


o o 




STATOR 


Figure 35. - LA survey locations used to investigate the 
periodically unsteady flow through a compressor 
stator. 




(a) Rotor trailing edge aligned circumferentially 
with stator leading edge. 



(b) Rotor displaced circumferentially by one-third 
rotor pitch. 



(c) Rotor displaced circumferentially by two-thirds 
rotor pitch. 


Figure 36. - Periodically unsteady stator flow field 
at three points in the rotor-stator blade passing 
cycle. Distribution of turbulence intensity 
parallel to local temporary main flow direction 
(Tpar) within one stator blade passage. 




2. Government Accession No. 


3. Recipient’s Catalog No. 


1. Report No- 

NASA TM-88798 


4. Title and Subtitle 

Laser Fringe Anemometry for Aero Engine Components 

5. Report Date 


6. Performing Organization Code 

505-62-21 

7. Author(s) 

Anthony J. Strazlsar 

8. Performing Organization Report No. 

E— 31 35 


10. Work Unit No. 

9. Performing Organization Name and Address 


' National Aeronautics and Space Administration 
Lewis Research Center / a//u 

Cleveland, Ohio 44135 ) 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Memorandum 

National Aeronautics and Space Administration 
Washington, O.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


Prepared for the 67th Symposium of the AGARD Propulsion and Energetics Panel on 
Advanced Instrumentation for Aero Engine Components, Philadelphia, Pennsylvania, 
May 19-23, 1986. 


16. Abstract 

Advances In flow measurement techniques In turbomachinery continue to be paced 
by the need to obtain detailed data for use In validating numerical predictions 
of the flowfleld and for use in the development of empirical models for those 
flow features which cannot be readily modelled numerically.) The use of laser 
anemometry In turbomachinery research has grown over the last 14 yr In response 
to these needs. Based on past applications and current developments, this paper 
reviews the key Issues which are Involved when considering the application of 
laser anemometry to the measurement of turbomachinery flowflelds. Aspects of 
laser fringe anemometer optical design which are applicable to turbomachinery 
research are briefly reviewed. Application problems which are common to both 
laser fringe anemometry (LFA) and laser transit anemometry (LTA) such as seed 
particle Injection, optical access to the flowfleld, and measurement of rotor 
rotational position are covered. The efficiency of various data acquisition 
schemes Is analyzed and Issues related to data Integrity and error estimation are 
addressed. Real-time data analysis techniques aimed at capturing flow physics 
In real time are discussed. Finally, data reduction and analysis techniques are 
discussed and Illustrated using examples taken from several LFA turbomachinery 
applications. 


17. Key Words (Suggested by Author(s)) 

Laser veloclmeter 

Turbomachinery 

Instrumentation 


19. Security Classif. (of this report) 

l 

20. Security Classif. (of this page) 

21 . No. of pages 

22. Price* 

Unclassified 

Unclassified 




*For sale by the National Technical Information Service, Springfield, Virginia 22161 


18. Distribution Statement 

Unclassified - unlimited 
STAR Category 02 

















National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland. Ohio 44 135 

Official Business 
Penalty for Private Use $300 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 


NASA 



